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A Robust Method for the Analysis of Integrated Spectra 
from Globular Clusters using Lick Indices 
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ABSTRACT 

We define a method for analysis of the integrated spectra of extra-galactic globular 
clusters that provides more reliable measures of age, metallicity and a-element abun- 
dance ratio than have so far been achieved. The method involves the simultaneous 
fitting of up to 25 Lick indices in a x^-fitting technique that maximises the use of the 
available data. Here we compare three sets of single stellar population (SSP) models 
of Lick indices to the high signal-to-noise, integrated spectra of 20 Galactic globu- 
lar clusters. The ages, [Fe/H] and a-element abundance ratios derived from the SSP 
models are compared to the results of resolved stellar population studies from the 
literature. We find good consistency with the published values, with an agreement of 
better than 0.1 dex in all three derived parameters. The technique allows identification 
of abundance ratio anomalies, such as the known nitrogen over-abundance in Galac- 
tic globular clusters, and the presence of anomalous horizontal-branch morphologies. 
It also minimises the impact on the derived parameters of imperfect calibration to 
the Lick system, and reduction errors in general. The method defined in this work is 
therefore robust with respect to many of the difficulties that plague the application of 
SSP models in general, and is, consequently, well suited to the study of extra-galactic 
globular cluster systems. 
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1 INTRODUCTION 

Even with the superior spatial resolution from space tele- 
scopes or adaptive optics, only for the very nearest galaxies 
will we be able to resolve individual stars. Thus the vast 
bulk of extra-galactic studies necessarily involve the analysis 
of integrated stellar populations. Models using Lick indices 
(Worthey 1994; Trager et al. 1998) have been developed 
which allow derivation of key parameters such as age, metal- 
licity and Qf-element abundance ratio from the spectra of in- 
tegrated stellar populations (e.g. Vazdekis 1999; Bruzual & 
Chariot 2003; Proctor & Sansom 2002; Thomas, Maraston & 
Bender 2003). These models are best suited to single stellar 
populations (SSPs), but they are often applied to the com- 
posite stellar populations found in galaxies (e.g. Worthey, 
Faber & Gonzalez 1992; Gorgas et al. 1997; Greggio 1997; 
Trager et al. 2000; Proctor & Sansom 2002; Terlevich & 
Forbes 2002; Proctor et al. 2004) 

Perhaps the best examples of single stellar populations 
that are easily accessible to us are the Galactic globular clus- 
ters (GCs), for which we have resolved colour- magnitude di- 



agrams (e.g. Piotto et al 2002) and/or high resolution spec- 
tra (e.g. Carney 1996). Indeed the SSP models are often 
built using such information from GCs. Thus, in order to 
have any confidence in the results of applying SSP models 
to unresolved galaxies, it is necessary that the SSPs correctly 
reproduce the properties of Galactic GCs. 

Here we apply three recent sets of SSP models to 
20 Galactic globular clusters with high signal-to-noise in- 
tegrated spectra. We also compare the results to those 
from high resolution spectra and colour-magnitude diagrams 
(CMDs) from the literature. Previous studies of Galactic 
GCs using Lick indices (Puzia et al. 2002; Thomas et al. 
2003) have used the GC ages measured from CMDs as an 
input to their interpretation of the indices. However, this 
approach is unsuitable for the studies of extra- galactic GCs, 
as such a priori knowledge of ages is unavailable. Instead, 
it is necessary to derive age, [Fe/H], and a-element abun- 
dance ratios directly from the spectra. In this way, we are 
able to test the SSP models and to place constraints on the 
accuracies that can currently be achieved for the ages and 
metallicities of extra-galactic GCs using Lick indices. 
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2 THE DATA AND MODELS 

We use data from two separate studies. The first is the data 
and spectra of 12 Galactic globular clusters (GCs) and the 
bulge of the Milky Way from the study of Puzia et al. (2002; 
hereafter P02). These long-slit data were collected using 
the ESO 1.52 m telescope on La Silla with a dispersion 
of 1.89 A/pix and achieved a signal-to-noise in excess of 
50 A~^. The observations covered a spectral range of ~3400 
to ~7300 A. This covers the Rose (1984) indices and 25 Lick 
indices from USa to the TiO indices (for Lick index defini- 
tions used see Trager et al. 1998 and Worthcy & Ottaviani 
1997). The main aim of the P02 study was the measurement 
of Lick indices in metal-rich bulge GCs. Eight such clusters 
were observed, as well as four associated with the halo. In 
order to ensure a representative sampling of the underlying 
stellar population a number of (generally three) individual 
spectra were obtained, offset from each other by a few arcsec- 
onds. The bulge spectrum was constructed from 15 different 
bulge fields, many in Baade's Window. Data reduction in- 
cluded calibration to the Lick system using 12 Lick standard 
stars observed on the same observing run. 

Our second source of Galactic GC data is from Cohen, 
Blakeslee & Ryzhov (1998; hereafter CBR98). Lick index 
values were taken from Beasley et al. (2004) . The data were 
obtained using the Keck I telescope with a dispersion of 
1.24 A/pix and achieved a signal-to-noise similar to that 
of the P02 data (i.e. S/N~50 A~^). For approximately half 
the GCs in the sample. Lick index measurements were made 
between 4200 A and 6500 A, while in the other half of the 
sample indices were measured between 4600 A and 6500 A. 
The study included twelve GCs from both the halo and 
bulge sub-populations. Table 1 shows that the sample con- 
tains red, blue and intermediate horizontal branch GCs. The 
combined P02 and CBR98 samples include 20 GCs (four of 
which were observed by both studies). 

Due to the lack of Lick standard star observations, the 
CBR98 data were not calibrated to the Lick system. For- 
tunately, there are four GCs in this data set common to 
the (fully calibrated) P02 study. The average differences be- 
tween the indices in the common GCs were therefore used 
to calibrate the CBR98 data to the Lick system. The off- 
sets were generally found to be ~0.3 A in line indices and 
~ 0.01 mag in molecular band indices, with no evidence of 
correlations with index values. These offsets are consistently 
smaller than the rms scatter in the differences between GCs. 
However, for Fe5015 an offset of 0.46 A (with 0.50 A rms) 
was found between the data sets (with the CBR98 data ex- 
hibiting lower values). We leave interpretation of this offset 
to Section 3.2. 

As this method of calibrating to the Lick system is far 
less reliable than using Lick standard stars, the index errors 
for the CBR98 data were taken to be the rms about the mean 
offsets for the four common GCs. This value is assumed for 
the error as it consistently exceeds the quoted observational 
error. However, we note that it is the relative size of index 
errors that influences the fitting procedure outlined in 
Section 3.2. 
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Table 1. Call index and horizontal branch ratios (HBRs; Harris 
1996) for Galactic globular clusters. 



2.1 Rose indices and horizontal branch ststrs 

We have measured the CaH index (Rose 1984) from the P02 
spectra. This index is known to be sensitive to the pres- 
ence of hot stars (Rose 1985) such as those occupying the 
blue horizontal branch (BHB). The measured CaH values 
are given in Table 1 and plotted against horizontal branch 
ratio (HBR) from Harris (1996) in Fig. 1. The three inter- 
mediate HBR GCs in this figure are shown by Piotto et al. 
(2002) to possess both red and blue horizontal branch stars. 
We shall therefore refer to such GCs as possessing interme- 
diate horizontal branches (IHBs). Fig. 1 suggests that the 
CaH index allows an estimation of the HBR in old stellar 
populations for which CMDs can not be obtained, i.e. from 
the integrated spectra of most extra-galactic systems. Schi- 
avon et al. (2004) detailed an alternative approach to the 
identification of GCs with BHBs, achieved by comparing 
the ratio of HSf to H/3 to their SSP models. We found the 
ratio predicted by the Schiavon et al. (2004) models to be 
less than the SSP model predictions used here (see Section 
2.2) for essentially all the GCs in our sample, i.e. the excess 
in the ratio in BHB GCs found by Schiavon et al. (2004) is 
not reproduced in the models used here. One possible ex- 
planation of this is that Schiavon et al. use the Jones (1996) 
stellar library, while the models used in this work are based 
on the Lick stellar library (see Section 2.2). 
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Figure 1. Call index against horizontal branch ratio (HBR) for 
the P02 sample. Intermediate horizontal branches (IHBs; trian- 
gles) arc evident in Piotto et al. (2002) colour— magnitude dia- 
grams of these GCs. Pure blue and red horizontal branches are 
shown by squares and circles respectively. The Call index pro- 
vides a good indicator of the HBR. 



2.2 The SSP models 

The results from three different SSP models are compared 
in this paper. These are: 

Vazdekis (1999; hereafter V99): The models used here 

arc those based on the fitting functions of Worthcy ct al. 
(1994) and Worthey & Ottaviani (1997). These SSP models 
cover the metallicity range -1.68 < [Fe/H] < 0.2 and ages 
from 0.1 to 17.78 Gyr. They use the isochrones of Bertelli et 
al. (1994). These isochrones include BHB stars in old, low- 
mctallicity populations. No account is taken of the variation 
in abundance ratios in the stars used to construct the SSPs. 
The index values therefore reflect the solar neighbourhood 
abundance ratio pattern (see Section 3.1), with a-element 
enhancements falling from ~-|-0.3 dex at low metallicities to 
0.0 dex at solar metallicity. 

Bruzual &c Chariot (2003; hereafter BC03): The mod- 
els used here are those derived from spectral energy distribu- 
tions and then calibrated to the lick system (i.e. from their 
lsindex_sedJick_system files). These SSP models cover the 
metallicity range -2.3 < [Fe/H] < 0.4 and ages from 0.1 Myr 
to 20 Gyr. The models use Padova (1994) isochrones (Alongi 
et al. 1993; Bressan et al. 1993; Fagotto ct al. 1994a,b; Gi- 
rardi et al. 1996). These isochrones include BHB stars in 
old, low-metallicity populations. No account is taken of the 
variation in abundance ratios in the stars used to construct 
the SSPs, so, again, the index values reflect the local abun- 
dance ratio pattern (see Section 3.1). 

Thomas, Maraston &c Bender (2003; hereafter 
TMB03): These SSP models cover the metallicity range 
-2.25 < [Fe/H] < 0.67 and ages from 1 to 15 Gyr. The 
models use the isochrones of Cassisi, Castellani & Castel- 
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Figure 2. The data of P02 and CBR98 are plotted on SSP grids. 
Symbols represent horizontal branch morphology as in Fig. 1. 

The average error in each index is indicated in the bottom right 
of each plot. Lines of constant metallicity range from [Fe/H]=0.0 
down to the lowest metallicity modelled (-1.5 for V99 and -2.25 
for BC03 and TMB03). Lines of constant age are from 1 Gyr up to 
the oldest population modelled (17, 20 and 15 Gyr for V99, BC03 
and TMB03 respectively). For TMB03 models, solar abundance 
ratio models are shown as solid lines, while models with an a- 
element enhancement of -1-0.3 dex are shown as dashed lines. The 
Galactic bulge (open circle) is also shown. The GCs NGC 6356, 
6528, 6553 and 6624 appear twice in each diagram as both the 
P02 and CBR98 values are presented. The data lie off the V99 
and BC03 grids at high index values, as they possess a-element 
enhancements greater than that modelled. 
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lani (1997), Bono et al. (1997) and Salasnich et al. (2000). 
TMB03 provides two versions of the Balmer lines in their 
SSP models; one modelling blue horizontal branches in low 
metallicity populations, the other modelling red horizontal 
branches (RHBs) at all metalUcities. TMB03 also adjust in- 
dex values to talce account of the variation in abundance 
ratios in the stars used to construct the SSPs using the re- 
suhs of Tripicco & Bell (1995; see Section 3.1). TMB03 thus 
provide index values for SSPs with both solar- and non- 
solar (a-element enhanced) abundance ratios. 

We refer the reader to the source papers for further de- 
tails. Here, we focus on obtaining reliable estimates of GC 
properties from each of the models. We note that Schiavon 
et al. (2002a,b) have produced SSP models, based on the 
colour-magnitude diagram of 47 Tuc. Unfortunately, these 
models are as yet unpublished and are, therefore, not con- 
sidered in this work. 

The data of P02 and CBR98 are plotted in the <Fe>- 
Mg2 plane in Fig. 2 where <Fc> = (Fe5270+Fc5335)/2. Agc- 
metallicity grids from the SSP models are also shown. Grid 
lines for age are shown from 1 Gyr up to the oldest age 
modelled. For metallicity, the grids show the lowest [Fe/H] 
modelled up to solar ([Fe/H]=0.0 dex). For TMB03 models, 
only the models with a -clement enhancement of +0.3 dex 
are shown. However, it is important to note that the TMB03 
models differ from those of V99 and BC03 in that TMB03 
correct their indices for the abundance ratio pattern in the 
stars used to construct the SSP models (sec Sections 3.1 
and 4.3). The TMB03 models, therefore, reflect a constant 
a-element enhancement at all metallicities, while V99 and 
BC03 SSP models reflect the varying abundance ratio with 
metallicity of the stars used in their construction. 

The effect of the varying abundance ratios with metal- 
licity in the V99 and BC03 models can be seen in Fig. 2. 
For all GCs, the TMB03 models suggest a near constant 
a-element enhancement which is consistent with the high 
resolution study of Carney (1996). However, in V99 and 
BC03 plots, while GCs with low <Fe> and Mg2 (i.e. low 
metallicity GCs), lie close to, or within, the grids, those 
with high index values lie well to the right of the grids. 
This is simply the result of low a-element enhancements 
in liigli-motallicity, local-neighbourhood stars used to con- 
struct the V99 and BC03 SSP models (e.g. Edvardsson et al. 
1993; Gustafsson et al. 1999; Bensby, Feltzing & Lundstr0m 
2003). 

The variations in abundance ratios severely complicate 
the process of deriving age and metallicity from the SSPs 
models. In order to overcome these difflculties, the [MgFe] in- 
dex was deflned as [MgFe]=^/ Mgb x (Fe5270 -|- Fe5335)/2 
(Gonzalez 1993). This combination of Lick indices was se- 
lected as it is insensitive to enhancement effects. Two com- 
monly used Balmer lines (H/3 and H7f) are plotted against 
[MgFe] in Fig. 3. Such two dimensional index plots are 
commonly used to derive ages and metallicities (e.g. Bres- 
san, Chiosi & Tantalo 1996; Kuntsclmcr & Davics 1998; 
Longhetti et al. 2000; Kuntschner et al. 2002; Beasley, Hoyle 
& Sharpies 2002). However, inspection of Fig. 3 clearly iden- 
tifles a number of concerns. Firstly, in both H/3 and U-fF 
plots, a significant number of the GCs fall below the oldest 
SSPs, suggesting ages significantly greater than the age of 



the universe. A second concern is that a number of GCs sug- 
gest young (^5 Gyr) ages, in conflict with the known age 
distribution of these Galactic GCs (Salaris & Weiss 2002; 
Rosenberg et al. 1999). In addition, for many GCs, the ages 
suggested by the H/3-[MgFe] grids differ from those sug- 
gested by the H7F-[MgFe] grids. For example, NGC 6171 
lies below the grids in H/3, suggesting an age >15 Gyr, while 
in the Hje plots its position within the grid suggests an 
age <5 Gyr. The bulge also suggests a younger age in the 
H7F-[MgFe] plots than in those of H/3-[MgFe] . On the other 
hand, the high metallicity GCs, NGC 6528 and 6553 (with 
[MgFe]~3.0 A), suggest older ages in the H7F-[MgFe] plots 
than in the H/3-[MgFe] plots. Therefore, the ages derived 
from Fig. 3 are sometimes highly dependent on the choice 
of Balmer line, and are inconsistent with the known age dis- 
tribution of these Galactic GCs, regardless of which Balmer 
line is used. Consequently, we conclude that 2-dimensional 
index plots, such as those in Fig. 3, are unreliable for the 
measurement of ages in globular clusters. In the following, 
we outline a method that is robust with respect to such 
problems. 



3 FITTING DATA TO SSP MODELS 

The method of measuring the derived parameters for GCs 

used in this work is a x^~niinimisation technique which in- 
volves fitting as many Lick indices as possible. The justifi- 
cations for this are: 

• By fitting as many indices as possible the procedure 
makes maximum use of the available data. It should be noted 
that, while individual indices are degenerate with respect 

to the derived parameters (age, metallicity and Q-element 
abundance ratios), most indices contain some information 
regarding each of these parameters. 

• The derived parameters are less prone to uncertainties 
in data reduction, e.g. flux calibration error, stray cosmic 
rays, skyline residuals, velocity dispersion error, errors in 
conversion to the Lick system, etc. 

• The procedure is also particularly useful for spectra af- 
fected by speciflc modelling problems, such as horizontal 
branch morphologies not matching those modelled in the 
SSPs. This is because the affected indices (mainly H/3 in the 
case of horizontal branch morphology problems) can be iden- 
tified and, if necessary, omitted from the fitting procedure 
with only a modest increase in overall error (see Proctor & 
Sansom 2002; hereafter PS02). Indeed, using this technique, 
the derived parameters are less prone to modelling uncer- 
tainties in general. 



3.1 Modelling the effects of non— solar abundance 
ratios 

In this work we use SSP models to derive estimates of 
age, [Fe/H] and the abundance ratio of a number of key, 
mainly a-, elements (parameterised by [E/Fe] - see PS02). 
TMB03 models include SSPs of solar abundance ratio (i.e. 
[E/Fe] = 0.0) as weU as [E/Fe] = -0.3, -1-0.3, -1-0.5. These 
models were constructed using the results of Tripicco & 
Bell (1995; hereafter TB95) by the method first outlined in 
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Figure 3. The data of P02 and CBR98 arc plotted on H/3-[MgFc] and H7i7-[MgFe] SSP grids. Symbols are as in Fig. 1. Grids lines in 
metallicity are from [Fe/H]=0.0 (solar) down to the lowest modelled (-1.5 for V99 and -2.25 for BC03 and TMB03) in steps of 0.25 dex. 
For age, grid Unas are 1, 1.5, 2, 3, 5, 8, 12 Gyr as well as the highest age modelled (17, 20 and 15 Gyr for V99, BC03 and TMB03 
respectively). There are fewer points in the H-fp plots, as this index is not measured in half of the CBR98 GOs. A wide range of ages is 
suggested by both H/3 and Hje plots, with many GCs falling outside the SSP grids. Large variations between the suggested ages in H/3 
and HjF plots are also evident, particularly for NGC 6171 and the bulge (open circle). 



Trager et al. (2000; hereafter TOO). V99 and BC03 SSP mod- 
els, on the other hand, have abundance ratios that vary with 
metalUcity, as their models reflect the abundance ratios of 
the stellar calibrators used to construct the models. We have 
calculated indices for V99 and BC03 SSP models for abun- 
dance ratios [E/Fe]=-0.3,-|-0.3 and -1-0.6, using the results 
of TB95. Two methods were considered. The first, (the TOO 
method), is best suited to low metallicity ([Fe/H]<0.0) pop- 
ulations. The second, described in PS02 (the 'Fe- method'), 
is better suited to high metallicity populations ([Fe/H]>0.0). 
Both of these methods are based on the assumption of a 
fixed fractional contribution to the total luminosity of each 
the three stellar types modelled by TB95. 

The fractional change in index values for a change of 
-1-0.3 dex in [E/Fe] axe listed in Table 2 for both methods. 
TB95 did not model H5 and H7 indices, so these indices axe 
assumed to be unaffected by changes in abundance ratios*. 



We note that, Thomas, Maraston & Korn (2004) have recently 



It is important to note that the two methods assume 
different chemical compositions, with the TOO method as- 
suming a fixed [Z/H] (overall metallicity) as [E/Fe] changes, 
while the PS02 method assumes fixed [Fe/H]. Unless other- 
wise stated, the derived parameters from these SSP models 
quoted in this work were derived using the TOO method. 
Linear interpolation (in grid steps of 0.025 dex) allowed us 
to estimate intermediate values of the derived parameters. 



3.2 Fitting procedure 

As our aim is to fit as many indices as possible, we be- 
gan by obtaining the best fit (by x^-minimisation) to all 25 



estimated the sensitivity of the H<5 and H7 indices to abundance 

ratios. However, they find that; . . the sensitivity of Hj and 
H6 to a/Fe disappears at low metallicity and is relatively small 
around solar metallicity ... ". 
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Figure 4. The average deviation in units of error (i.e. x) for GCs in the P02 sample. Indices are listed in wavelength order. Error bars 
represent the rms scatter in the deviations. Open circles show the results of including all 25 indices in the fitting procedure. Solid circles 
represent the deviations when poorly fitting indices (open squares) are omitted from the fitting procedure (see text). 



observed Lick indices in log(agc), [Fe/H] and [E/Fc] param- 
eter space. For the TMB03 SSP models, fits were obtained 
to both the RHB and BHB models. 

Residuals to the best fit (observed value minus best fit 
value expressed in terms of index errors, i.e. x) ^rc sum- 
marised in Fig. 4. This figure shows the average and rms 
deviations from the best fit values for all indices of the P02 
sample, and for all three sets of SSP models. It is clear from 
this figure that a number of indices are poorly fit by all 
three sets of SSP models. The poorest fits were obtained for 
the CN indices and NaD. The NaD index is known to be 
severely affected by interstellar absorption. This index was 
therefore excluded from the fitting procedure. On the other 
hand, the CN indices, which were also found to be enhanced 
by Trager (2004), are extremely sensitive to nitrogen abun- 
dance (TB95) which is known to be higher in Galactic GCs 
than in field and halo stars (e.g. Li & Burstein 2003). The 
large positive residuals present in the fits to GCs (but not 
evident in the fit to the P02 bulge data) is therefore con- 
sistent with an [N/Fe] > in GCs. Indeed, TMB03 showed 
that SSP models in which N is independently enhanced with 
respect to other light elements cam provide excellent fits to 
these indices. The fitting procedure has therefore success- 
fully identified this abundance anomaly. Rather than intro- 



duce N abundance as another free parameter in our mod- 
elling, the CN indices were excluded from the fitting proce- 
dure in the following. 

As the CN and NaD indices were removed form the fit- 
ting procedure, the fits of other indices (particularly the en- 
hancement sensitive indices C24668, Mgi, Mg2 and Mgb) 
are generally improved (see Fig. 4). However for two indices 
(Fe5Q15 and Ca4227) this was not the case. For Fe5015 we 
recall that this index showed a large variation between the 
two samples for the four GCs in common to P02 and CBR98. 
The positive deviation in Fig. 4 is consistent with the as- 
sumption that P02 overestimate the Fe5015 index values. 
In this case the fitting procedure therefore appears to have 
identified a problem in calibration to the Lick system. Con- 
sequently, Fe5015 was excluded from the fitting procedure. 

The Ca4227 index was also shown by TB95 to be sensi- 
tive to nitrogen abundance (with increased N abundance re- 
sulting in a reduction in the Ca4227 line— strength) . Its neg- 
ative deviation in Fig. 4 is therefore consistent with the ni- 
trogen overabundance in GCs. However, Ca4455 also shows 
a negative deviation. This weakness in the calcium lines has 
previously been observed in GCs (e.g. Maraston et al. 2003) 
and galaxies (e.g. Worthey 1992; Vazdekis et al. 1997; Trager 
et al. 1998; PS02), as well as in the near- infrared CaT in- 
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Index 


PS02 


TOO 




method 


method 


HSa 






HSf 






CNl 


0.058 


0.029 


CN2 


0.071 


0.032 


Oa42z7 


-0.092 


-0.272 


G4300 


0.097 


0.053 


H7A 






H7F 






TTl O C5 O 

J?e4383 


-0.029 


-0.132 


Ca4455 


0.185 


0.012 


Fe4531 


0.054 


-0.077 


Fe4668 


0.368 


0.032 


Hbeta 


0.031 


0.031 


Fe5015 


0.059 


-0.077 


Mgl 


0.048 


0.024 


Mg2 


0.048 


0.018 


Mg b 


0.315 


0.207 


Fe5270 


-0.012 


-0.127 


Fe5335 


-0.083 


-0.196 


Fe5406 


-0.058 


-0.172 


Fe5709 


0.015 


-0.068 


Fe5782 


0.049 


-0.159 


Na D 


0.258 


0.025 


TiOl 


0.002 


0.001 


Ti02 


l).0()2 


-0.001 



Table 2. Fractional change in indices for an [E/Fe] change of 
-1-0.3 dex, calculated from TB95 (see PS02). TB95 did not model 
H5 and H7 indices. These indices are therefore assumed to have 
no sensitivity to [E/Fe]. It should be noted that application of 
the two methods differs as they describe different chemical com- 
positions (see PS02). 

dex (e.g. Saglia et al. 2002). It is therefore unlikely to be a 
result of problems in calibration to the Lick system. In addi- 
tion, since TB95 showed Ca4455 to be insensitive to calcium 
abundance, it is unlikely that this is the effect of a variation 
in [Ca/Fe]. Finally, we note that the indices either side of 
Ca4455, particularly Fe4383 and Fe4531, also possess con- 
sistently negative deviations, suggesting that there is some 
systematic discrepancy between models and observations at 
these wavelengths. Wo therefore conclude that the poor fit of 
Ca4455 is most likely the result of problems in the modelling 
of this index in SSP models. 

Despite its relatively large deviations in Fig. 4, exclusion 
of Ca4455 had negligible impact on the results of the fitting 
procedure, with 85% of the fits finding exactly the same 
best fit, while the remaining 15% differed by at most 2 grid 
steps (amounting to <0.05 dex in any derived parameter). 
To be consistent with our aim of fitting as many indices as 
possible wo therefore decided to include this index in the 
fitting procedure. 

Having considered the average and rms deviations to 
the best fits, we turn our attention to the deviations of in- 
dices in individual GCs. For the most part, no single index 
dominated the deviations in any given GC. The main ex- 
ception to this being H/3 in seven GCs fitted using either 
V99 or BC03 SSP models, as well as the TMB03 RHB mod- 
els. For the three of these aberrant GCs in the P02 data. 
Call values are all below 1.15, indicating that these are all 



either BHB or IHB GCs. The decision was therefore taJsen 
to exclude H/3 from all GCs in the P02 sample with Call 
below this value, i.e. in all the BHB and IHB GCs. For the 
four remaining CBR98 GCs with deviations dominated by 
H/3 deviations, the index was also omitted from the fitting 
procedure. It is worth noting that these four GCs are also 
those in the CBR98 data with HBR values >0.9 and there- 
fore constitute all the BHB and IHB GCs in the CBR98 
sample. 

The decision to exclude certain H/3 values is supported 

by the finding of PS02 that the exclusion of Balmer lines 
from the fitting procedure has relatively little impact on the 
derived parameters, as long as a sufficiently large number 
of indices are included in the fitting procedure. The results 
from TMB03 SSP models are taken to be those from the 
RHB fits in GCs with HBR=-1.0, and for BHB fits for GCs 
with HBR>-1.0 (i.e. IHB and BHB GCs). 

Finally, the deviations from the best fits of TiO indices 
in NGC 6626 and G4300 in NGC 6637 and NGC 6981 were 
also found to be large. Consequently, these indices were ex- 
cluded from the fitting of the respective GCs. The final de- 
viations from the best fits to the P02 data are shown in Fig. 
4. 

A similar process was carried out for the CBR98 data. 
However, these data spanned a narrower wavelength range 
and, consequently, it was only possible to include between 8 
and 16 indices in the fits (compared to between 18 and 20 
using P02 data) . For common indices the results were similar 
to those for the P02 data, with the exception of the TiO 
indices which were found to give poor fits in the CBR98 data. 
TiOi and Ti02 were therefore excluded from the fitting of 
the CBR98 GCs. H7 indices in NGC 6171 were also found 
to possess large residuals to the best fits (see Fig. 3) . These 
indices were therefore omitted from the fitting of this GC. 

The best fits of the P02 data to all three SSP model 
sets exhibited average reduced-x^ values of ~3.5. This sug- 
gests that either reduction errors exist that have not been 
accounted for, and /or that significant errors are present in 
the SSP models. Indeed, the typical reduced values for 
our fits to CBR98 data (for which the rms scatter in the 
comparison of 4 GCs common to both P02 and CBR98 was 
used as the error) were ~1.5, suggesting that both effects 
noted above are present. However, whatever the case, our 
use of a large number of indices in the fitting procedure 
ameliorates the effects of such problems on the estimates of 
derived parameters. It is also worth noting that, since the 
reduced of the best fits were significantly greater than 
1.0, results for lower signal-to-noise spectra can be expected 
to be similar. However, from previous experience (e.g. Proc- 
tor et al. 2004), we find that a minimum signal-to- noise of 20 

is required for the technique to obtain robust results. 

In order to test the robustness of the fitting procedure, 
fits were obtained for each of the GCs with only eight indices 
included in the fitting procedure (we recall that between 
eight and 20 indices were included in the final fits above, 
with the majority of GCs fit by 16 or more indices). The 
eight indices included in these fits were H/3, Mgi, Mg2, Mgb, 
Fe5270, Fe5335 Fe5406 and Fe5709. The rms difference be- 
tween the derived parameters from these fits and those from 
the final fits above were 0.103 dex, 0.144 dex and 0.087 dex 
in log(age), [Fe/H] and [E/Fe] respectively. 

In summary, the x^-minimisation technique has allowed 
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us to identify problems with NaD (probably caused by inter- 
stellar absorption) , the known nitrogen abundance anomaly 
affecting CNi, CN2 and (possibly) Ca4227, problems in cali- 
bration to the Lick system (Fe5015 in the P02 data and TiO 
indices in CBR98) , as well as individual indices sufTering un- 
explained reduction errors. This technique therefore holds a 
great deal of promise in the investigation of integrated stellar 
populations such as extra-galactic globular cluster systems. 



4 COMPARING DERIVED PARAMETERS TO 
LITERATURE VALUES 

In this section we review the derived parameters (age, 
[Fc/H], [Z/H] and [E/Fe]) obtained from the fitting proce- 
dures described in the previous section. These arc denoted 
with the subscript SSP. The results are compared to those 
obtained by other techniques (such as fitting to the H/?- 
[MgFe] and H7F-[MgFe] grids shown in Fig. 3) and to values 
from the literature. We wish to emphasise that the values 
of the derived parameters were not considered in the fitting 
procedure outlined in the previous section, i.e. the 'best' fits 
were arrived at without considering the derived parameters 
that they generate. 

Errors in the derived parameters were estimated us- 
ing Monte-Carlo type realisations of the best fit SSP mod- 
els perturbed by the index error estimates. However, this 
method produces errors in the derived parameters based 
purely on the estimates of observational error, they do not 
include any allowance for modelling errors. We estimate 
more realistic errors in Section 5. 

Six GCs (NGC 6121, 6205, 6218, 6341, 6626 and 6981) 
were found to possess [Fe/H] ssp below the minimum mod- 
elled by V99. These GCs have therefore been omitted from 
presentations of V99 results. NGC 6341 also fell below the 
minimum [Fe/H]ssp in BC03 and TMB03 SSP models. This 
GC has therefore been omitted from presentations of BC03 
and TMB03 results. 



4.1 Age estimates 

Before discussing the results of the age determinations as a 
whole, it is first necessary to consider the method employed 
in the fitting of the high metallicity GCs NGC 6528 and 
NGC 6553, as these possess mctallicitics close to solar, i.e. in 
the range affected by the changing sensitivity of isochrones 
to non-solar abundance ratios (see PS02). It was shown in 
PS02 and Thomas & Maraston (2003) that the main effect of 
this on the values of derived parameters is in age. We confirm 
this finding using the GC data of both P02 and CBR98, as 
values of [Fe/H] ssp and [E/Fe] ssp are consistent between 
the TOO and PS02 methods''' (within errors) for both GCs in 

t Wc note that the application of the PS02 method to TMB03 
solar abundance ratio SSP models causes a small (unavoidable) 
internal inconsistency in these models, as the TMB03 models are 
already corrected for the local abundance ratio pattern ■using a 
method similar to TOO. However, the effects of this are expected 
to be negligible, as the local abundance ratio, for which TMB03 
compensate in generating their solar abundance ratio SSPs, is 
close to the solar value for SSPs with [Z/H]~0.0 dex. 



both studies. Values of age, on the other hand, show signifi- 
cant differences. The average age and rms value (from both 
GCs compared to three sets of SSP models) using the PS02 
method is 10.3±1.5 Gyr, while values derived by the TOO 
method were 7.8±3.3 Gyr with one value as low as 4 Gyr. 
For the bulge, the PS02 method finds ages of 6.0±0.8 Gyr 
while the TOO method gives 3.8±0.4 Gyr. The PS02 method 
therefore produces age estimates that are higher than those 
derived by the TOO method (as noted in PS02 and Thomas 
& Maraston 2003), and more consistent with the known old- 
ages of GCs (Safaris & Weiss 2002; Rosenberg et al. 1999) 
and the bulge (e.g. Rich 1999). Wc have therefore elected 
to use the derived values obtained by the PS02 method for 
NGC 6528, 6553 and the bulge in the following analysis. 

For the four GCs common to the P02 and CBR98 stud- 
ies (NGC 6356, 6528, 6553 and 6624) we find an average 
age offset of only 0.02 Gyr and an rms scatter of 1.7 Gyr 
about a one-to-one relation, although we note that the two 
studies arc not fully independent, as the CBR98 data were 
calibrated to the Lick system using a comparison to the P02 
data (Section 2). 

Ages measured from the SSP models for 10 GCs 
are compared with published ages derived from colour- 
magnitude diagrams (CMDs) in Fig. 5. CMD-bascd values 
are taken from Safaris & Weiss (2002) for NGC 6121, 6171, 
6205, 6218, 6356, 6624, 6637 and 6838. For NGC 6528 and 
6553 we use the Bruzual ct al. (1997) result that these GCs 
are ~2 Gyr younger than 47 Tuc, which is given as 10.7 Gyr 
old in Salaris & Weiss. We therefore assume an age of 8.7 Gyr 
for NGC 6528 and NGC 6553. 

Fig. 5 shows GC ages derived in three ways. The top 
plots show ages derived by interpolation of the H/3-[MgFe] 
grids (Fig. 3) . The tendency for GCs to fall below the grids 
in Fig. 3 is reflected by arrows indicating ages greater than 
or equal to that modelled. The extremely young ages (3- 
5 Gyr) suggested for some GCs (primarily IHB and BHB 
GCs) is also evident. Indeed, there is generally poor agree- 
ment with ages derived from CMDs. The middle plots show 
ages estimated by interpolation of the H7F-[MgFe] grids 
(Fig. 3). Again, a large scatter in ages, and poor agreement 
with ages derived from CMDs are found. The inconsistency 
of age estimates derived from comparison to the H/3-[MgFe] 
and H7F-[MgFe] grids, noted in Section 2.2, is evident, par- 
ticularly for NGC6171, as well as NGC 6528 and 6553 (both 
with agecMD <10 Gyr). This emphasises the need for a 
method of age measurement that is robust to both observar 
tional errors and modelling deficiencies. 

In the bottom plots the results of our minimisation 
technique are shown. In this case, no GCs are found to pos- 
sess ages less than 8 Gyr, including IHB and BHB GCs 
which exhibit old ages. In addition, and despite the narrow 
age range suggested by CMD values, we find correlations be- 
tween log(age)ssp and log(agc)cMD at 95%, 95% and 90% 
confidence levels for V99, BC03 and TMB03 SSP models 
respectively (the correlation in TMB03 results are adversely 
affected by the relatively low upper age limit of the models). 
The x'^-minimisation technique therefore produces age esti- 
mates that are in reasonable agreement with values derived 
from CMDs in both absolute and relative terms. We there- 
fore conclude that this technique is superior to the method 
of deriving ages from 2-dimensional index plots such as Fig. 
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Figure 5. Ages derived from H/3-[MgFc], H7i7-[MgFc] and the x'^~rninimisation technique (top, middle and bottom plots respectively) 
are compared to values derived from CMDs. Values of a,gccMD from Salaris & Weiss (2002) and Bruzual ct al. (1997). One-to-one lines 
are shown in each plot (solid lines), as well as the highest SSP age modelled (dashed line). Points that lie off the Balmcr-[MgFe] grids are 
plotted at the highest age modelled and identified with an arrow. GCs for which the x^^minimisation technique found ages equal to the 
oldest modelled are also identified by an arrow. While very poor agreement is found with CMD ages for the Balmer-[MgFe] grid methods, 
relatively good agreement is found using the x^-minimisation technique. Note also the age variations between the Balmer-[MgFe] grid 
methods. 



3 and is, consequently, well suited to the estimation of ages 
in extra-galactic GC systems. 

4.2 Metallicity estimates 

In this section we compare the metallicities derived from 
SSP models ([Fe/H]ssp and [Z/H]ssp) with values pub- 
lished in the literature. 

For the four GCs common to the P02 and CBR98 stud- 
ies (NGC 6356, 6528, 6553 and 6624) we find an average 



[Fe/H] difference of only 0.01 dex and an rms scatter of 
0.08 dcx about a one to one relation, although wc again 
note that the two studies are not fully independent, as we 
have calibrated the CBR98 data to the Lick system using a 
comparison to the P02 data. Given the good agreement be- 
tween studies average values are presented in the following. 

In their analysis of P02 data, TMB03 found a good 
correlation between [Z/H] and the [Fe/H] values of Harris 
(1996). A comparison of our [Fe/H]ssp and [Z/H]ssp mear 
sures, from V99, BC03 and TMB03 SSP models, to the val- 
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Figure 6. Metallicitics ([Fc/H] and [Z/H]) derived from SSP models are compared to the [Fc/H] values given in Harris (1996). Symbols 
shapes represent horizontal branch morphology as in Fig. 1. Symbol size represents GCs from P02 (large symbols) and CBR98 (small 
symbols). A good agreement is achieved between the SSP derived value of [ZfKjssp and [Fe/H] from Harris (1996) for all three sets of 
SSP models. 



ues of [Fe/H] of Harris (1996) are shown in Fig. 6. The same 

good agreement between [Z/H] ssp and the [Fe/H] values 
of Harris (1996) is found for all three sets of SSP models, 
with rms variations about a one-to-one lines <0.2 dex. How- 
ever, we note that Beasley et al. (2002) find good agreement 
between [Fe/H] derived from SSP models and values from 
CMDs for GCs in the Large Magellanic Cloud. We therefore 
simply conclude that the methods used in this work (and 
TMB03) are well suited to the study of metallicities in extra- 
galactic GC systems, as they produce excellent measures of 
relative metallicities. They also successfully recover the ac- 
tual [Fo/H] value assigned to the isochrono that best fits the 
observed population - i/ wo assume that [Z/¥L]ssp rcfiects 
the [Fe/H] of the isochrone. However, this result suggests 
that there are still some problems in the calibration of ei- 
ther the P02 indices (which we used to calibrate the CBR98 
indices to the Lick system) or the metallicity scale in SSP 
models. 



4.3 Abundance ratio estimates 

In this section we compare the values of [E/Fe]ssp, derived 
by the minimisation technique, with values of enhance- 
ments in Q-elements published in the literature. 

It is important to boar in mind that, for V99 and BC03 
SSP models, the values of [E/Fe] output from the fitting pro- 
cedure are enhancements relative to that of the best fit SSP 
models (see PS02, TMB03 and Proctor et al. 2004 for more 



complete discussions). Consequently, since the SSP models 

are constructed using stars in the solar neighbourhood, ac- 
count must be taken of the pattern of enhancements in local 
stars (see Section 3.1). TMB03 have adjusted their indices 
for the [E/Fe] inherent to local stars such that [E/Fe] is 
quoted with respect to solar at all metallicities. The TMB03 
SSP models therefore require no further correction. How- 
ever, for V99 and BC03 the output values of [E/Fe]gsp need 
to be corrected for this effect. In this work we assume a pat- 
tern in local abundance ratios (and hence in the SSP models) 
similar to that assumed by TMB03. The [E/Fe] of the SSP 
models is assumed to rise from 0.0 dex at [Fe/H]ssp=0.0 dex 
to [E/Fe]=+0.25 dex at [Fo/H]s5p= 1.0 dex and to remain 
at this values for all [Fc/H]ssp < 1.0 dex (e.g. Edvards- 
son et al. 1993; Gustafsson et al. 1999; Bensby, Feltzing & 
Lundstr0m 2003). 

A comparison of a-enhancements from the literature 
to [E/Fe]ssp for five GCs (NGC 6121, 6205 and 6838 from 
Carney 1996, NGC 6528 from Carretta et al. 2001 and 
NGC 6553 from Cohen et al. 1999) is shown in Fig. 7. For 
the sake of consistency with Carney (1996), we use the aver- 
age of the Si, Ca and Ti abundances to calculate a-element 
euhauccnicnts from the latter two studies. The agreement 
between [E/Fe] ssp and the published value of [Si,Ca,Ti/Fe] 
for these five GCs is good. We therefore conclude that the 
X^-minimisation technique provides an excellent technique 
for estimating abundance ratios in extra-galactic GC sys- 
tems. 
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[E/Fe]Lit 

Figure 7. A comparison of [E/Fc] derived from SSP models to 
those from the literature. Symbols shapes represent HBR as in 
Fig. 1. Symbol size represents GCs from P02 (large symbols) and 
CBR98 (small symbols) . The one-to-one lines are also shown. The 
agreement between the two measurements is good. 



5 DISCUSSION 

Given that our aim is to define a procedure for measuring the 
derived parameters of extrargaJactic GCs, it is clearly useful 
to ask: What conclusions would we have drawn from the data 

as a whole had it been an extra-galactic GC system? Can 
this technique help us in the investigation of such systems? 
To this end, the results of our determinations of the derived 
parameters for the 19 CCs studied are plotted in Fig. 8. 

First, we note that the fitting procedure outlined in Sec- 
tion 3.2 identified the known abundance anomaly affecting 
CN indices (assumed to be the result of nitrogen abundance 
enhancement). The capability of the technique to highlight 
such abundance variations could clearly provide useful in- 
sights into extra-galactic CC systems. 

With regard to age estimates, it was shown that the 
X^-minimisation technique provides values that are in good 
agreement with CMD derived ages in both an absolute and 
relative sense. The average values of age (± the rms scat- 
ter) were found to be 13.1±2.3, 12.2±3.3 and 12.7±1.9 Cyr 
for V99, BC03 and TMB03 respectively. Therefore, all three 
SSP models used indicate an old population of GCs. Given 
that at least some of the scatter about these averages is the 
result of real variations in GC age, the rms scatters about 



the average values of age indicate that we can safely assume 
an upper limit of 0.1 dex on errors in individual log(age) 
estimates. Such a value is also consistent with the rms de- 
viations found in the test of robustness detailed in Section 
3.2. 

The metallicity determinations for all three SSP models 
show a broad range of sub-solar metallicities, which are in 
good agreement with the literature values (see Section 4.2). 
The scatter about the one-to-one relations in the top plots 
of Fig. 6 indicate that we can safely assume an upper limit of 
±0.1 dex on errors in individual [Fe/H] and [Z/H] estimates. 
This value is again consistent with the rms deviations found 
in the test of robustness detailed in Section 3.2. 

The [E/Fe] SSP estimates for the five GCs with lit- 
erature values for a-enhancements are in good agreement 
with the literature values (Section 4.3), again in both an 
absolute and relative sense. The average and rms values 
of [E/Fe]ssp for the sample are 0.35±0.09, 0.25±0.10 and 
0.32±0.09 dex for V99, BC03 and TMB03 respectively. We 
can again use the rms values to set an upper limit on the er- 
rors in [E/Fe]ssp of ±0.1 dex, which is again consistent with 
the rms deviations found in the test of robustness detailed 
in Section 3.2. 

Interestingly, Fig. 8 suggests that values of [E/Fe] ssp 
in GCs are slightly higher than that of the local stars used 
to populate the SSPs at all metallicities. The additional 
enhancement in a-elements, relative to that of the stellar 
calibrators, is 0.14, 0.05 and 0.11 dex in V99, BC03 and 
TMB03 SSP models respectively. It should be noted that 
these [E/Fe] excesses are independent of the assumed [E/Fe] 
pattern in the stellar calibrators, as the fitting technique out- 
puts the [E/Fe] relative to the models. While an enhance- 
ment with respect to the stellar calibrators is perhaps to be 
expected at high metallicities ([Fe/H] > -1), where the SSP 
models are based largely on thin disc stars, the effect persists 
to the low metallicities at which the models are constructed 
largely using thick disc and halo stars. This is a slightly 
surprising result that clearly requires confirmation. 

The stellar bulge data appear to indicate a relatively 
young age (6.3±3.6). However, the P02 bulge data are highly 
susceptible to stochastic fluctuations in the stellar contents 
of the fields from which the spectra are constructed. Given 
this, and the large error on the age estimate, no conclusion 
is drawn here regarding the age of the bulge. On the other 
hand, the average [Fe/H]ssp of -0.37 dex and [Z/H]ssp of 
-0.09 dex found for the bulge are consistent with the known 
distribution of bulge star metallicities (Ferreras, Wyse & Silk 
2003), while the [E/Fe]ssp of 0.29 dex is consistent with the 
a element enhancements from high resolution studies (e.g. 
Mc William & Rich 2004). Therefore, while no firm conclu- 
sion is reached with respect to age, the x^~niinimisation 
technique produces metallicity and abundance ratio esti- 
mates for the bulge which are in good agreement with liter- 
ature values. 

In summary, we find that the values of age, [Fe/H] 
and [E/Fe] derived using the x'^ - minimisation technique are 
therefore in good agreement with values published in the 
literature and reproduce many of the known properties of 
Galactic GCs and the bulge. Using the data of P02 and 
CBR98, upper limits on statistical errors in log(age), [Fe/H] 
and [E/Fe] have been shown to be ±0.1 dex. Derived param- 
eters with such errors could clearly provide powerful insights 



© 0000 RAS, MNRAS 000, 1-13 



12 Proctor, Forbes & Beasley 

V99 



BC03 



TMB03 



D_ 
(/) 

to 

CD 

o 



(/) 
C/) 



O 



O 



LO 



0.6 
17 0.4 

Li_ 

\ 
LjJ 

^ 0.2 





T — I — I — I — I — I — r 



1 TTT 1 ItT 1 1 1 1 1 1 1 1 n 1 1 1 1 1 1 1 1 1 1 1 r 




-2 -1 
[Fe/H]3sp 



-2 -1 
[Fe/H]33p 



-2 -1 
[Fe/H]ssp 



Figure 8. The values of age and local abundance ratio corrected [E/Fe] are plotted against [Fc/H] derived from SSP models. Symbol 
shapes represent horizontal branch morphology as in Fig. 1. Symbol size represents GCs from P02 (large symbols) and CBR98 (small 
symbols). Values for the bulge are also shown (open circles). There are no discernable differences between the results for the two datasets. 
Plots of age against [Fe/H] (upper) include the average age of the sample as horizontal lines. The x^^mhiimisation technique finds GCs 
to be consistently old. The [E/Fe] plots (lower) show the assumed local abundance ratio pattern as solid lines. GCs appear slightly 
enhanced in a-elements compared to the SSP models at all metallicities. The bulge exhibits a similar [E/Fe] to the GCs. 



into the chemical and formation histories of extra-galactic 
GC systems and their host galaxies. 



rameters. We therefore conclude that the x -minimisation 
technique holds groat promise for the study of the ages and 
metallicities of extra-galactic GC systems. 



6 CONCLUSIONS 

We have compared the Lick index data of Puzia et al. (2002) 
and Cohen ct al. (1998) for 20 Galactic GCs to the SSP 
models of Vazdekis (1999), Bruzual & Chariot (2003) and 
Thomas, Maraston & Bender (2003). We showed that ages 
derived from 2-dimcnsional index plots, such as H/3-[MgFe], 
are not reliable for these data. An alternative approach of 
using a x^-minimisation technique was outlined. This tech- 
nique allowed the identification of the known nitrogen over- 
abundance in Galactic GCs, as well as interstellar absorp- 
tion in NaD and poorly calibrated indices. It was shown 
that this technique successfully recovers the known age, 
mctallicity and abundance ratio properties of the Galac- 
tic GC population, reproducing values from the literature 
with an accuracy of ~0.1 dex in any parameter. However, 
even after poorly modelled indices were omitted from the 
^^-minimisation procedure, the fits to the models still ex- 
hibited reduced-x^ values significantly in excess of 1.0 for 
all GCs in both studies. This implies that SSP modelling 
uncertainties axe sigiuficant, and emphasises the advantages 
of using a large number of indices in measuring derived pa- 
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